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Abstract 



We investigate the properties of single-particle resonances in a non- 
spherical potential by solving the coupled-channels equations for the radial 
wave functions. We first generalize the box discretization method for positive 
energy states to a deformed system. As in the spherical case, we find that 
the discretized energy is stabilized against the box size when a resonance con- 
dition is met. Using the wave functions thus obtained, we then discuss the 
energy and the radial dependences of scattering wave functions in the vicinity 
of an isolated resonance. In the eigenchannel basis, where the 5-matrix is 
diagonal, we propose a generalized expression for the factorization formula 
for the multi-channel wave function. We find that the factorized wave func- 
tion agrees well with the exact solution inside the centrifugal barrier when 
the energy distance from the resonance is less than the resonance width. 



1 



Typeset using REVTfrX. 



I. INTRODUCTION 



Recent experimental activities on drip-line nuclei have opened up a renewed interest in 
the continuum spectra of finite many-body systems. These nuclei are often weakly bound, 
or may even be unbound, and a proper treatment of the positive-energy states is crucial 
for their theoretical description. The continuum influences the properties of loosely bound 
nuclei mainly through the resonant states, which have a dominant contribution both in the 
ground state, e.g., through the correlations induced by the virtual pair scattering, and in the 
excited states [1-9]. It is thus becoming more and more important to approach the structure 
of unstable nuclei from the nuclear reaction point of view [10,11]. Recent developments of 
shell model calculations for weakly-bound systems by Id Betan et al. [1] and by Michel et 
al. [2] clearly indicate this direction. 

A single-particle resonance state can be obtained relatively easily when the mean field 
potential has a spherical shape. In this case, the single-particle angular momentum is 
conserved, and the Schrodinger equation is reduced to a one dimensional differential equation 
for the radial motion. Approximated resonance states are easily found by putting a system 
in a box and imposing a condition that the wave function vanishes at the edge of the box. 
With this boundary condition, the continuum spectra are discretized. It has been recognized 
that, for a resonance state, the single-particle energy is stabilized against the box size, while 
it changes significantly as a function of the box radius for non-resonant continuum states 
[12-14]. 

A more consistent treatment for the resonance state is also possible for a spherical system. 
The resonance states are associated with poles of the S'-matrix in the complex energy plane. 
One way to obtain the resonance states is to impose the outgoing wave boundary condition 
in the asymptotic region with a complex energy. The state which satisfies this boundary 
condition is referred to as the Gamow state in the literature. A similar approach is used in 
the complex scaling method, which has been developped originally in atomic and molecular 
physics [15] and subsequently been applied to nuclear systems [16,17]. Alternatively, the 
resonance states can also be found by monitoring the energy derivative of the phase shift 
along the real energy axis [18], where a maximum appears at the resonance energy. For a 
scattering wave function in a spherical potential, a simple factorization formula has been 
derived for its energy and radial dependences in the vicinity of an isolated resonance state 
[18-20]. This factorization property has been used recently in order to estimate the effect 
of resonant states on pairing correlations in loosely bound nuclei [5]. 

A resonance structure is much more complicated for a deformed system, where several 
angular momentum components of the single-particle wave function are coupled to each 
other. A standard method to solve this problem is to expand the wave function on a discrete 
basis, which may be deformed harmonic oscillator wave functions [21,22] or eigenfunctions 
of a spherical potential [23,24]. The positive energy states obtained by diagonalizing the 
Hamiltonian in such basis can describe properly the main properties of the resonant states 
[13,14], but this may require a large number of the basis set. A more direct method to 
obtain resonance states in a deformed system is to setup the coupled-channels equations 
for the radial motion and solve them for a complex energy by imposing the outgoing wave 
boundary condition at infinity for all the open channels [25]. Recently, this method has 
been successfully applied to the proton emission decay of proton- rich nuclei [26,27]. This 
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method has an advantage over the matrix diagonalization method that the wave functions 
are subject only to the radial mesh discretization error for the Schrodinger equation and 
that the desired asymptotic boundary condition can be imposed independently of the size 
of the mesh interval. However, the description of resonant states with such complex wave 
functions is generally difficult to use in nuclear models, since all the evaluated observables 
become complex quantities if the rest of the complex non-resonant states is not taken into 
account. We thus prefer to keep the energy as a real variable here. 

The aim of this paper is to develop a practical method to obtain scattering wave functions 
in the vicinity of a multi-channel resonance state. A difficulty is that there are N linearly 
independent solutions of the coupled-channels equations at a given energy E (N is the 
number of included channels), whose asymptotic behavior are all different. For example, a 
physical scattering boundary condition is defined by requiring that asymptotically there is 
an incoming wave only in the incident channel. In this case, the equations have N degenerate 
solutions depending on which channel is the incident one. Since any linear combination of 
these N solutions is also a solution of the coupled-channels equations, a problem arises as to 
which combination one should take in order to study efficiently a multi-channel resonance 
state. If the resonance width is extremely small, as in the case of proton decay, one can 
easily construct a resonance wave function by matching to a standing wave solution [28-30]. 
However, this procedure is not applicable to a broad resonance. 

Our method presented in this paper is to start with discretized states whose wave func- 
tions vanish at a given radius. This is an extension of the box boundary condition for a 
spherical system into the deformed regime. As in the spherical case, resonances can be 
identified with those states whose energy is stabilized against the box size. This method 
provides the most convenient basis states, where the unperturbed states are already a good 
approximation to the true resonance states. In fact, as we will show below, the resonance 
energy thus estimated is close to the energy at which the first derivative of the eigenphase 
sum with respect to energy has a maximum. Here, the eigenphase sum is defined as the 
sum of the phase shifts for the eigenchannels, for which the S- matrix is diagonal [31]. In 
the field of electron-molecule scattering, it has been known that the eigenphase sum has 
the same energy dependence as the elastic phase shift in a one-channel case [32,33]. Using 
these unperturbed wave functions, we then employ the first order perturbation theory to 
construct wave functions around a resonance state. 

The paper is organized as follows. In Sec. II, we discuss how to generalize the box 
boundary condition to a multi-channel case. This uses a generalized concept of node for 
a multi-channel wave function [34]. We also discuss the method based on the eigenphase 
sum and compare it with the multi-channel box discretization method. In Sec. Ill, we 
introduce the eigenchannels and discuss their usefulness. We consider several cases: the 
uncoupled (spherical) case, a very narrow resonance, and a case where the coupled-channels 
equations are decoupled. In Sec. IV, we discuss the factorization of the wave functions in 
the vicinity of an isolated resonance. We first consider a perturbation to the box discretized 
wave functions for a spherical system in order to re-formulate the factorization formula of 
Unger [19], which was originally derived using the i?-matrix theory [10,35-37]. We then 
apply it to the eigenbasis wave functions, and derive a generalization of the factorization 
formula for a multi-channel wave function. We examine the validity of the factorization 
formula by comparing it with the exact solutions of the coupled-channels equations. The 



3 



summary of the paper is then given in Sec. V. In the Appendix, we present an explicit 
form of the coupled-channels equations for a Hamiltonian which contains a central and a 
deformed spin-orbit potentials. 



II. SINGLE-PARTICLE RESONANCES AND BOX-DISCRETIZED STATES 

Let us start with the following single-particle (s.p.) Hamiltonian: 

ff= _*W(r), (1) 

where the potential V may also contain a derivative operator such as a spin-orbit interaction. 
When the potential V is non-spherical, the s.p. angular momentum is not conserved and 
the total s.p. wave function has the form 

*(r)=E^^>W*)> (2) 

Ijm 

with 

yijrn{r) = { lm i\ m s\jm)Y lmi (r)Xm s - (3) 

mi,m 3 

Here, Y\ mi are spherical harmonics and \m s denotes the spin wave function. Projecting the 
Schrodinger equation, H^/ = E^/, on the spin-angular states the coupled-channels 

equations for the radial wave functions then read, 



h 2 d 2 l(l + l)k 2 „ 
— E 



2m dr 2 2mr 2 



u ljm (r) + r J2 (yijmlVly^m')^^ = 0. (4) 

I'j'm' r 



These equations are solved with certain boundary conditions which depend on the prob- 
lem of interest. A standard way to solve them is to generate N linearly independent solutions, 
N being the number of channel states to be included, and then to take a linear combination 
of these N solutions so that the relevant boundary condition is satisfied [28,38]. The lin- 
early independent solutions can be obtained by taking iV different sets of initial conditions 
at r = 0. We denote these solutions by <pLL'{f), where L refers to the channels while V refers 
to a particular choice of initial conditions. Here, we use a shorthand notation, L = (Ijm). 
A simple choice for the N initial conditions is to impose 

<f) LL ,(r)^r l+1 5 LyL , for r -> 0. (5) 

The coupled-channels equations are solved outwards for each initial condition V . The wave 
functions Ul(t) in Eq.(4) are written in terms of the <Pll' as 

UL(r) = £CW LL ,(r), (6) 
v 

where the coefficients Cy are determined from the asymptotic behavior of u L (r). 
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When the s.p. energy E is positive, the physical wave function is a scattering wave with 
an incoming wave behaviour in some particular incident channel L and outgoing waves in 
all channels L. Thus, the asymptotic boundary condition of the wave function ul{t) is given 
by 



Mr) - \ { e ~ t(kr ~ ln/2)5 L,L - W i(fcr ~ W2) } for r -> oo, (7) 



where k = y2mE/h 2 and Sll is the scattering S'-matrix. The normalization is chosen 
so that the total wave function satisfies {^e,l \^e,l' ) — — E')S LojL > q . Here, |^_b,l ) * s 
given by Eq. (2), whose channel components Ul(t) satisfy the boundary condition (7). In 
practice, the S-matrix can be obtained by decomposing <pLL'(f) in the asymptotic region 
with spherical Hankel functions h^ftr) and comparing Eq. (6) with Eq. (7) [28,38]. 

The generalization of the concept of s.p. resonances from the spherical to the multi- 
channel cases is most conveniently done by monitoring the energy dependence of the eigen- 
phase sum [32,33]. The eigenphases are related to the eigenvalues of the S'-matrix through 

(WSUU, = e 2i5 °8 a , a , (8) 

Hazi [33] has shown that, the sum of the eigenphases A(E) = J2 a $a(E), has the same energy 
dependence around a resonance as the phase shift in a spherical system , i.e., 

A(E) = A (E) + tan- 1 2{E ^_ E y (9) 

where the slowly- varying quantity A (E) is the sum of the background eigenphases, and E R 
and T are the multi-channel resonance energy and total width, respectively. Therefore, the 
properties of multi-channel resonances can be extracted directly by plotting A(E), or its 
energy derivative, as a function of E. The quantity A(E) is called the eigenphase sum, and 
it has been widely used in the context of electron- molecule scattering (see, e.g., Ref. [39] for 
a recent publication). 

We illustrate the above discussion by the following example. We use a Woods-Saxon 
parametrization for the potential V which is given by 

V(r) = Kcnt(r) + V(V, s (r)) • (-iV x <x), (10) 

with 

KentW = V (r) - R p^l Y20 (r), (11) 

dV (r) 

VUr) = V^{r)-R so (3 a ^^Y 2Q {r), (12) 

where V (r) and V so (r) have a Woods-Saxon shape: 

V (r) = -V /[l + exp((r - R)/a)}, (13) 
V so (r) = V so /[l + exp((r - R so )/a so )]. (14) 
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For simplicity, we have expanded the deformed Woods-Saxon potential and kept only the 
linear order of the deformation parameter (3 [40] . Also, we assume an axial symmetry, where 
both the parity it and the spin projection K onto the z axis are conserved. An explicit 
form of the coupled-channels equations is given in the Appendix. The parameters of the 
Woods-Saxon potential are taken to be Vo=42.0 MeV, R = R so = 1.27 x44 1/3 fm, a = a so 
= 0.67 fm, and V^, = 14.9 MeV/fm 2 , to somehow simulate the neutron potential in the 44 S 
region. As an example we calculate the neutron s.p. levels with K n = 5/2 + at (3—0.2. We 
include the d 5 / 2 , #7/2, #9/2 ,213/2, and iu/ 2 states in the coupled-channels equations. Notice 
that the si/ 2 and rf 3 / 2 states do not contribute to the K 7 " = 5/2 + levels. We have checked 
that the results are not significantly altered even if we include higher angular momentum 
components. 

The thick solid line in Fig. 1 shows the eigenphase sum and its energy derivative as a 
function of energy. The contributions from three main eigen-channels are also shown by thin 
lines. One clearly sees a maximum at E R — 3.44 MeV in the first derivative of the eigenphase 
sum. The total width can be estimated to be T= 0.46 MeV. 

Let us now compare this continuum result with the s.p. spectrum of a box-discretized 
calculation. From the structure (6) of the general solution, it is clear that one can impose 
the boundary condition that the total wave function \&(r) vanishes at a radius -Rbox- This 
is equivalent to putting the nucleus in an impenetrable spherical box of radius -Rbox- If the 
condition 



is satisfied, then Eq.(6) has a non-trivial solution for Cl such that ul(R^, ox ) = for all 
the channels L. This is a natural extension of the well-known box discretization method 
for a spherical system. We notice that Johnson [34] has advocated to use det(0(r)) as a 
generalized concept of node for a multi-channel wave function. We have recently used this 
method to find all the bound state solutions of a deformed Skyrme-Hartree-Fock mean field. 

Figure 2 shows s.p. energies obtained with the box discretization method as a function of 
the box size -Rbox- One can clearly see that there are two classes of s.p. levels. One consists 
of those whose energy changes significantly as the box radius -Rb ox increases, and the other 
contains those whose energy is almost constant as a function of -Rbox- This behavior is 
similar to that of box-discretized s.p. energies for a spherical system where the stabilized 
levels correspond to resonances [12]. This example shows that, in the deformed case it is also 
possible to generate a complete set of box-discretized states and that some of these states 
can be a good approximation to multi-channel s.p. resonances. 



The states obtained with the box discretization method in the previous section are the 
eigenstates of a potential which is the same as the original potential for r < -Rbox, and 
infinite for r > -Rbox- They thus form a complete set, and any regular function defined in 
the domain r < .Rbox can be expanded on this basis. In the next section, we discuss the 
factorization property of the scattering wave function around a resonance in terms of this 
basis set. To this end, it appears convenient to introduce the eigenchannel wave functions 
[31] defined in terms of the unitary matrix U of Eq. (8), 




(15) 



III. EIGENCHANNEL REPRESENTATION 
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(16) 



We refer to the channels L as the physical channels, in order to distinguish them from the 
eigenchannels a. Substituting Eqs. (7) and (8), one can find that the eigenchannel wave 
functions behave asymptotically as 

{e-^-eV-^)}^ jm (f) for r ^ oo. (17) 

For each eigenchannel, the asymptotic radial wave functions thus behave in the same way 
for all angular momentum components, and the coupling is greatly simplified. Loomba et 
al. argued that the eigenchannel wave functions provide the most direct way to visualize 
a resonance state [31]. The eigenchannels were also used in the context of the R matrix 
theory [37,41,42]. We also mention that the eigenchannel approach has been widely used in 
the field of heavy- ion fusion reactions in the context of fusion barrier distributions [43-47]. 
In this section, we discuss the usefulness of this approach for multi-channel resonances by 
considering several cases. 

A. Spherical systems 

In spherical systems, the s.p. angular momentum is conserved and the S- matrix is 
diagonal in the physical channel representation L, i.e., the physical channels are equivalent 
to the eigenchannels. Likewise, the wave function (piuij') with the initial condition (5) is 
diagonal at all the values of r, and Cl in Eq. (6) is proportional to 8l,l for the incident 
channel L . From this consideration, it is apparent that the eigenchannel wave functions ty a 
exhibits a resonance behavior for a spherical system if the resonance condition is met for an 
incident channel a = L . 

B. Two-channel decoupled system 

We next consider a special case where the coupled-channels equations are completely 
decoupled by a transformation of the basis. For this purpose, we consider the following 
two-channel model [43], 

(4£-<HCT + U) T)G£j)-* ™ 

This corresponds to a system where the incident channel couples to a spinless vibrational 
state whose excitation energy is zero [43,44,47]. This model has been used by Dasso et 
al. in order to understand the coupling-assisted tunneling phenomena in heavy-ion fusion 
reactions [43]. 

We refer to the basis in the coupled-channels equations (18) as |1) and |2). For symmetry 
reasons the S- matrix in this system satisfies Su = S 2 2 and S±2 = S 2 i, and the total wave 
functions are given by 
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*i = ^Mr)|l)+« 2 (r)|2)], (19) 

* 2 = -Mr)|l)+«i(r)|2)]. (20) 
r 

Here, the asymptotic behaviors of the wave functions U\{r) and u 2 (r) are 

ui(r) ^ e- lkr - S n e ikr , (21) 

and 

u 2 (r)^-S 12 e ikr , (22) 
respectively. By diagonalizing the S matrix, the eigenchannel wave functions read 

*i = + *2) = l -Mr) + u 2 {r)) ■ (23) 

*2 = " *2) = ;(«i(r) - «2(r)) • 1^^. (24) 

Notice that the coupling matrix in the coupled-channels equations (18) can be diagonalized 
independently of r with the same unitary operator U. Transforming Eqs. (18) with this 
unitary operator leads to 

Ln ' " 2 +VIt)-e)(^ t \ +v ^\) + (^ °, M U fl +U2( , r l)=0. (25) 



2m dr 2 J \ui(r) — u 2 {r) J \ -F(r) J \ui(r) - u 2 {r) 



The coupled-channels equations are thus decoupled with this basis. The radial wave func- 
tions for the eigenchannel wave functions \^i j2 obeys a one-dimensional Schrodinger equation 
with a potential given by V(r) ± F(r), respectively. Therefore, if the one-dimensional po- 
tential V(r) +F(r) holds a resonance state at an energy En, the eigenchannel wave function 
exhibits a resonance behavior while the wave function ^2 does not, and vice versa. 
The energy dependence of the wave functions with the physical basis ^2 is much more 
complicated. Although, in general cases, the coupled-channels equations are not decoupled 
if different angular momentum components are coupled to each other, it is evident from 
this example that the eigenchannel basis provides a much simpler image of a multi-channel 
resonance than the physical basis. 



C. Very narrow resonance 



As a final example, we consider a case where the resonance width is very small, as in 
proton or alpha decays. In Refs. [28-30], the energy of a multi-channel resonance for a proton 
emission decay was found by imposing the asymptotic boundary condition that the channel 
wave functions are proportional to Gi(kr) for all the channels, where Gi(kr) is the irregular 
Coulomb wave function. From Eq. (17), one can see that this state should correspond to 
one of the eigenbasis wave functions whose eigenphase is n/2. In this subsection, we would 
like to demonstrate that the total width is exhausted by only one eigenchannel when it is 
very small. 
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To this end, let us start with the Breit-Wigner expression for a mult i- channel ^-matrix 
[10,35-37]. If there is only one isolated resonance of energy Er in the energy region of 
interest, the S'-matrix can be aproximated in the vicinity of this resonance as 



S w = e^5 L , L , - e^ + M (26) 

where 0l is the background phase shift and T is the total width given as a sum of partial 
widths T^. For a very narrow resonance, the background phase shifts can often be neglected. 
The left hand side of Eq. (8) then reads 

where = Y.l V^lU i a is related to the partial width for the eigenchannel a. By definition 
of the unitary operator U, the off-diagonal components of the right hand side of Eq. (27) 
should vanish. The only possibility for this to be satisfied is that T a is zero except for 
one particular channel a . Thus, around the resonance the eigenphase shifts are given by 
8 ao = tan _1 [r/(2(_E^ — E))] for the eigenchannel a , and 5 a = for the other channels. 
Again, this example clearly shows that the eigenchannel representation provides the most 
direct manifestation of a multi-channel resonance state. 



IV. FACTORIZATION OF SCATTERING WAVE FUNCTIONS 

We now consider the most general case for a multi-channel resonance, i.e., when the 
resonance width is not very small. For a spherical system having an isolated single-particle 
resonance at energy Er, Unger [19] has shown that the scattering radial wave function at 
energy E whose asymptotic form is 



ip E (r) -> J -^-e i5{E) sin(A;r - ln/2 + 5(E)) for r -> oo, (28) 
V tiE 

can be approximately factorized into a product of an r-independent function F(E) times 
the wave function at resonance energy. This factorization property holds if 

\E-E R \ « \V(v)\, (29) 

a condition which is realized generally in the internal region of the potential when \E — Er\ 
is less than the resonance width V. In Ref. [19] the function F(E) is written in terms of a 
Breit-Wigner function, but for an isolated resonance it is easy to relate this to the derivative 
of the phase shift with respect to energy, so that the factorization property can be expressed 
as 



f R dr\i> ER (r)\i 



under the condition (29). Thus, the factorization function F(E) is just e lS<yE ^ ^ jr| except for 
a normalization factor. In this section, we generalize this formula to a deformed system. We 
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do not attempt to give an analytic derivation, which could be done by using multi-channel 
Jost functions [48], but we would like to have a hint on how should be the generalized form 
of Eq.(30). To this end we find it useful to use the box discretized states of Sec. II to express 
the function F(E), first in the familiar spherical case, then in the general case. 



A. Spherical system 

The infinite set of all discretized states X x (r) with corresponding eigen-energies E x cal- 
culated in a box of radius -Rbox form a complete orthonormal set, and one can expand on this 
basis any regular function defined in the domain [0,i?box[ (the domain excludes r = i?box)- 
In particular, the inner part of the scattering wave function ipE{ r ) can be expanded as 

Mr) = T, C x(E)X x (r) re[0,R hox [, (31) 

A 

with 

fRbox 

C X (E) = drX x (r)^ E (r). (32) 

Jo 

Notice that this expansion does not hold at r = R hox in general. Here and in the following, 
all integrals J flbox dr are meant to be lim e ^o Jo Rb °*~ e dr. With this definition, the expansion 
(31) is still meaningful even when the radius r is infinitesimally smaller than i?box- Since 
the wave functions X x {r) and ipE( r ) obey the same Schrodinger equation (but with different 
energies from each other), 

n*^ E + ( v{r) + i(ip)V_ E]i , E = 0t 



2m dr 2 

Ti 2 d 2 X x ( 1(1 + l)Ti 



2m dr 2 

the expansion coefficients C\(E) are given by 




+ V(r) + v ' ' -E x \X x = 0. (33) 



Ca(£) = — ft oZT / dr ^e-tt - x. 



E-E x 2mh V dr ' 2 dr 2 J : 

= E~^ x 2^ MRhox)X ' x{Rhox) ■ (34) 

In the last step we have used the boundary condition X x (R hox ) = 0. 

In the vicinity of Er ~ E Xr , the contribution from A = \r dominates the expansion (31). 
Retaining only this term in the expansion (this corresponds to the one term approximation 
in the .R-matrix theory [10,35-37]), the scattering wave function around the resonance and 
inside the potential barrier reads 

Mr) - Y^ET ' ^M^ ox )X' XR (R hox )X XR (r). (35) 
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We now set the i?-matrix radius R to be -Rbox, approximate X\ R (r) by the resonance wave 
function ipE R { r ) (except for the normalization) and finally assume that /c-R box is sufficiently 
large so that the ^(r) at _Rb ox can be replaced by its asymptotic form (28). Then, a 
comparison between Eqs. (30) and (35) leads to 



e* 5iE) \^ = T^V • elHE) Sin ^* " + ^)) r X ' XR(Rh ° X) TV2 



E — E 



(36) 

We have thus established the relation between the function F(E) (except for a normalization 
factor) and the derivative of the discretized state X\ R at r = Rbox- This relation is our 
guideline for obtaining the generalized function F(E) in the deformed case. Notice that we 
have explicitly written down the normalization of X Xr (r) so that Eq. (36) can be used also 
when the function X\ R {j-) is not normalized. 

B. Deformed system 

The procedure is very similar to that of the spherical case. Using the complete set of 
discretized wave functions introduced in Sec. II : 

*A(r)=£^W), (37) 
we expand the eigenchannel wave functions ^E, a (f)' 

L ' 

= £C A (£)X A (r), (38) 

A 

with 



= £/ °*drX LX (r)ut\r) . (39) 



a. 

L 



The coupled-channels equations (4) lead to 

Cx(E) = • |^E4 a) (^ox)^ A (i? box ) . (40) 

As in the spherical system, we retain only the contribution from A = \r in the expansion 
(38) for E ~ E\ R . In the asymptotic region (for a large -Rbox) the wave function (r) has 
exactly the same asymptotic form (except for a factor Ulo) as the scattering wave function 
-0e(r) in the spherical system (see Eqs. (17) and (28)), the only difference being that the 
spherical phase shift 5(E) is replaced by the eigenphase 5 a . Thus, we can replace in Eq.(40) 
■u^(-Rbox) by ipE(Rbox)ULa with ipE represented by its asymptotic form (28). Comparing 
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the resulting expression with Eq. (36), we obtain the desired factorization formula for a 
multi-channel resonance, 



i6a(E) t 



ndE 



£ 



dr\X LXR (r) 



1/2 



u 



La 



(41) 



r 1 1/2 

Here, J2l $o h °* dr X\ R {r) 2 U La is the probability of finding the particle inside the box 
for the eigenchannel a. 

In practice, one may replace the box wave function X\ R in Eq. (41) by the scattering 
wave functions at the resonance, for each eigenchannel a, i.e., 



*E, a (r) ~ C a (E)* ER , a (r)/C a (E R ) 



(42) 



with 



Ca(E) = indE-- 



1/2 



U La {E) 



El fo Rb °* dr \u[ a \r)\i 



1/2 



(43) 



where the wave functions u 



(a). 



are evaluated at the resonance energy One may also 



choose the phase factor of the eigenchannel wave functions ^E,a f° r each channel a and 
energy E so that the factor e t5a does not appear in Eq. (43). With this, the wave functions 
become real numbers. 

Figure 3 shows the applicability of the factorization formula for the system considered 
in Figs. 1 and 2. At the resonance energy, E = Er = 3.44 MeV, the eigenphase shifts 
5 a /ir are calculated as 0.559, 2.53xl0" 3 , 4.12xl0~ 2 , 1.06xl0~ 5 , and 2.79xlCT 5 for the five 
eigenchannels (see the thin lines in Fig. 1). Therefore, there is one eigenchannel which 
dominates the eigenphase sum A(E) around the resonance. In the figure, the d 5 / 2 , #7/2, 
and gg/2 wave function components are shown for this particular eigenchannel. The exact 
wave functions at the resonance energy are denoted by the thick solid line, while those 
at E = 3.54 MeV are given by the thin solid line. The dashed line is the result of the 
factorization formula (42) and (43). We take i?b ox = 20 fm to evaluate the integrations in 
the formula. For the / = 4 components, one clearly sees that the energy dependence, as well 
as the radial dependence for r < 10 fm, of the radial wave functions are well expressed by 
the factorization formula. We have confirmed that the agreement is even better for the / = 6 
components (not shown). For the / = 2 component, the wave function is not localized inside 
the barrier, and the factorization formula does not work well. As expected, the factorization 
approximation is better for high I states, which are more localized by the centrifugal barrier. 

Figure 4 shows the spin projected total wave function (multiplied by r), 



(Xi/2\rV E ,a{r)) = £ 



u 



Ij 



[r)-(l,K-± i, ±\jK)Y lK _ 1/2 (r), 



(44) 



for this eigenchannel a. The upper and the lower panes are the wave functions as a function 
of r for {Q,4>) = (7r/4,0) and (Q,4>) = (n/2,0), respectively. (Note that the total wave 
function vanishes along 9 = for K = 5/2.) As we have shown in fig. 3, the factorization 
formula works reasonably well, especially inside the potential barrier, i.e., r < 4.23 fm. 
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In Fig. 5, we show more in detail the energy dependence of the main component of 
the multi-channel wave function (i.e., the g 9 / 2 component) around the resonance energy, 
and its comparison with the factorization method. At E — E R + 0.1T, the factorization 
formula works, and the exact radial dependence of the wave function is well reproduced, as 
we already showed in Fig. 3. As the energy increases, the performance of the factorization 
becomes worse. At E = Er + T, the factorization still works inside the potential, but the 
deviation becomes significant in the outside (see the middle panel). By the energy Er + 2r, 
the wave function does not resemble the resonance wave function, and the factorization 
formula looses its applicability, as was already known for the spherical system [20]. 

In Ref. [49], Lejeune and Mahaux used a solvable two-level model to argue that the energy 
dependence of the wave function components is quite different from each other and a simple 
scaling does not hold. However, this does not come out in the present study. The energy 
dependences are similar for all the angular momentum components, except for the <i-wave 
which shows a non-resonant behavior. This different conclusion may come from the fact 
that here the calculations are done with the eigenchannel basis, in contrast to Lejeune and 
Mahaux, who worked with the physical basis. This fact again suggests that the eigenchannel 
basis provides a powerful means to study a multi-channel resonance. 

V. SUMMARY 

We have used the eigenchannel representation of the multichannel wave function to in- 
vestigate the properties of a single-particle resonance in a deformed system. We first showed 
that the box boundary condition can be generalized to deformed systems if one uses the 
generalized concept of the wave function node, that is to impose the condition that the 
determinant of the wave function matrix vanishes at the boundary of a spherical box. We 
have solved the coupled-channels equations with this boundary condition and showed that 
the discretized energy is nearly a constant as a function of the box size for a multichannel 
resonance state. The energy E R and the width V of the resonant state were also calculated 
by using the energy dependence of the eigenphase sum. We have next studied the prop- 
erties of the eigenbasis which diagonalizes the S'-matrix. We argued that the eigenchannel 
representation has the most direct connection to a multi-channel resonance state. Using 
this representation, we have derived a formula which relates the wave functions at different 
energies around a resonance state, and studied its applicability. For an isolated resonance, 
we have observed a nice agreement between the approximate and exact wave functions inside 
the potential radius at energies \E — Er\ < T. 

In many-body problems, one often takes the single-particle states as a good building block 
to start with. This is the case, for instance, in the BCS calculation for a pairing Hamiltonian. 
If one includes unbound states in the calculation, one then has to deal with integrals over 
a wide range of energies in order to compute matrix elements. In the vicinity of a narrow 
resonance, the energy interval has to be taken small, which may be quite time consuming. 
This difficulty is substantially lessened if one uses the factorization formula derived in this 
article. Work is now in progress which uses this formula for resonance contributions in 
a Hartree-Fock plus BCS calculation for deformed nuclei. We will report the results in a 
separate paper. 
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APPENDIX A: COUPLED-CHANNELS EQUATIONS WITH A DEFORMED 

SPIN-ORBIT POTENTIAL 



In this Appendix, we present an explicit form of the coupling matrix elements in the 
coupled-channels equations, Eq. (4), for a Hamiltonian with a deformed spin-orbit interac- 
tion in the form of Eq.(lO). We first perform the multipole decomposition for the potentials: 

V CCIlt (r) = J2V c ^\r)Y x ,(r), (Al) 
W = £^(r)*V(r). (A2) 

We use the Wigner-Eckart theorem [50] to evaluate the matrix elements. For the central 
potential U cen t, this yields [51] 

/-*> IT> |-V> •,u r i' m '{.r) 
f\yijm I ^ccnt | yj'l'm' ) " 

= E(-r m (j m I ^)(^dln||^)U c ^(r)^ w (r), (A3) 



where / is defined by y2l + 1. 

For the spin-orbit part, we first decompose it as [50] 



V(U ls (r)) • HV x a) = £ { -\l ^ " ^\r) ] [*W-<V x <r^ 



+ \/2ATT (^P + ^^"Mj ln-(-iV x „)]<*> } , (A5) 
and use the Wigner-Eckart theorem for the operators [Y x ±i(—iV x <x)] A ^. That is, 

.Uj'l'm'ir) 



rOhjmWV^r)) - (-W x <r)\y. 



j'l'm' 



E(-r m 



-m /i m' 
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X < -1 



A + l dV, 



(A M ) 



A, 



- y^\r) ) r(y sl \\\Y x+1 (-iV x <r)]W\\y fl ) 



2A + 1 \ dr 



Uj'l'm'{r) 



+1 



A / dV, 



A + l, 



+ — ^is V) M ) r(^||[n-i(-^V x <r)]W\\y fl ) 



2A + 1 \ dr 



Uj>i'm'{r) 



The reduced matrix elements are given by [51] 

(y fl \\\y^i(-iv x tr)ny jV 
= E(-) AW+1 v^w{ A ^ 

r z" r n 

1/2 1/2 1 



Ail 1 A 



x < 



I / / 1 J 



J 3 

(^■ill^ll^»i»><Hlv||0, 



with 



(/"|| v| |0 



VF+T(|:-f) if l" = V + l, 

-y/v(i + ^) if i" = i'-i- 



This yields 



r(^||[lA ± iHVx CT )]( A )||^) 



Uj'l'm'(r) 



l'+3/2 



r 

Ail 1 A 

3 f 



j"=l> +1/2 L ^ 

' + 1 V 1] , , 

1/2 1/2 1 (^||1a||^ +1 ) • y/F+1 



d r + r 



Ufl'm'(r) 



3 

V-3/2 



' 1 



i E (-) AW+I v / 36W{ A ^ 

[V-l V 1 

1/2 1/2 1 

3" 3' 1 



Ail 1 A 

3 3" 



x < 



(^HIaII^-i)-(-v^) + 



Uj'l'm'(r) 



(A6) 



(A7) 



(A8) 



(A9) 
for the 



The matrix elements of the coupling potential were evaluated in this way in Ref. 
vibrational excitation in the proton decay of spherical nuclei. 

Notice that the matrix elements become trivial for the monopole part (A = 0), and are 
given by 



r(yi jm \Vx =0 \y 



j'l'm' 



Uj'l'm'(r) 



-(00), 



v£2(r) + I d\ I ■(/■) 



An 



Anr dr 



(A10) 
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1.5 




E (MeV) 

FIG. 1. The eigenphase sum (the upper panel) and its energy derivative (the lower panel) as a 
function of energy for neutron positive energy states of 44 S at (5 = 0.2. The spin projection onto 
the symmetry axis and the parity are taken to be K n = 5/2 + . These are denoted by the thick solid 
line, while the contributions from three main eigenchannesl are also shown by the thin dashed, 
dotted and dot-dashed lines. 
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FIG. 2. Neutron single-particle energies of 44 S obtained with the box discretization method as 
a function of the box radius -Rbox- The axial symmetry is assumed, with the deformation parameter 
[3 of 0.2. The spin projection onto the symmetry axis and the parity are taken to be K n = 5/2 + . 
The potential parameters are the same as in fig. 1. 
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FIG. 3. The channel wave functions around the neutron resonance state of S at (3 = 0.2 
and for K w = 5/2 + . Those are for the eigenchannel basis which dominates the total width. The 
thick and thin solid lines are the wave functions at the resonance energy, Er = 3.44 MeV and at 
£7=3.54 MeV, respectively. The dashed line is obtained by using the factorization formula for a 
multi-channel resonance state. The potential parameters are the same as in fig. 1. 
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FIG. 4. The total eigenchannel wave function around the neutron resonance state of 44 S at = 
0.2 and for K w = 5/2 + . The spin degree of freedom is eliminated by projecting the wave function 
onto the upper spin component. The upper and the lower panels are for the wave function along 
the (9,(f>) = {tt/4, 0) and the (6, (ft) = (tt/2,0) directions, respectively. The meaning of each line is 
the same as in fig. 3. 
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FIG. 5. The energy dependence of the #9/2 component of the eigenchannel wave function around 
the K w = 5/2 + resonance of 44 S. The solid line is the exact wave function, while the dashed line 
is obtained with the factorization formula for a multi-channel wave function. 
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